Characterisation of high energy electron irradiation damage in UPt 3 samples 
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We present transport and specific heat measurements on high quality single crystals of UPt3 
before and after irradiation by high energy electrons. We observe a strong dependence of the critical 
temperature with the sample thickness. The dramatic effects of the irradiation on the specific heat in 
the superconducting state can then be simply explained in terms of an inhomogeneous distribution 
of superconducting transition temperatures. The question of the failure of the "universal limit" 
in the heat transport of UPt3 is reexamined, and the conditions for a clean experimental test are 
established. 



PACS numbers: 

I. INTRODUCTION 

Unconventional superconductors are mainly found 
among strongly correlated electronic systems, where im- 
purities seem to be dominantly in the unitary limili. In 
such a case, even a very small concentration of defects will 
induce a finite density of states at zero energy^. More- 
over, a finite impurity concentration induces conduction 
bands which contribute to transport properties^: for ex- 
ample, impurities should add a linear term in tempera- 
ture (T) to the thermal conductivity for T— > 0. More 
subtle effects were predicted depending on the particular 
gap topology. If a line of nodes is present on the Fermi 
surface, and in presence of defects, the linear term of the 
thermal conductivity should be universal ( i.e. indepen- 
dent of the impurity concentration) for a small variation 
of T,4£. This is considered as a very robust test of the 
unconventional nature of the superconducting state and 
of its gap topology. 

Such a universal limit has been observed in a 
few quasi-2D unconventional superconductors such 
as YBa 2 Cu 3 6 .g&, Bi 2 Sr 2 CaCu 2 0gI and Sr 2 RuOA 
But it failed in the underdoped High-T c compound 
YBa 2 Cu40gS. In this last case, the origin of the ab- 
sence of the universal limit is still not clear, and could 
be attributed among others to the nature of the super- 
conductivity or the particular ground state of this com- 
pound. Further investigations have to be done in com- 
pounds where the electronic ground state is well estab- 
lished. Moreover, the universal limit has never been ob- 
served in any three dimensional unconventional super- 
conductor. 

One of the best studied unconventional 3D supercon- 
ductor is the heavy fermion UPt3. Its superconducting 
state is particularly interesting because it exhibits a well 
known phase diagrami°*ii, each phase corresponding to a 
different symmetry of the superconducting gap. Various 
scenarios exist for these phases^, corresponding to differ- 
ent types of nodes of the superconducting gap. Probing 
the universal limit in the different crystallographic direc- 



tions would be a very strong test of the models^. 

Impurity effects in UPt3 have been studied in many 
different ways. Substitution of Uranium or Platinum by 
another element have been used to study the decrease 
of the critical temperaturaiiiii^, and the interplay be- 
tween the magnetism and the superconductivity^. They 
would not be appropriate to probe the universal limit due 
to the rapid change of the normal state properties and T c 
with even the smallest concentrations that can be reason- 
ably studied. Small defect concentrations can be tuned 
with annealing, which varies the number of intrinsic met- 
allurgical defects. A main advantage of annealing is the 
possibility to follow the effects on the same sample. Its 
influence on physical properties at low temperature like 
the specific healiii or transport properties^ has been re- 
ported. Nevertheless, it does not seem appropriate to 
investigate a possible universal limit in UPt3 due to the 
small range of change of the critical temperature and the 
lack of knowledge of the nature of the defects : extended 
defects are expected rather than point defects, and they 
may require more involved descriptions 18 which spoil the 
strength of the theoretical predictions. 

A third technique uses high energy electron irradiation 
damage to create point defects^. It seems ideal for such 
studies, combining the advantages of the previous tech- 
niques : the vacancy and interstitial defects created by 
electron irradiation have a size of a few Angstroms, small 
compared to the coherence length of UPt3 (£o ~ 100 A); 
Moreover, it is possible to access a wide range of defect 
concentrations simply by tuning the irradiation time. Fi- 
nally, the same sample can be studied for the various 
concentrations, and a simple annealing at 700° C sup- 
presses again all defects, refreshing the sample. This 
technique has been used in UPt3 to study the influence 
of impurities on T c , on the upper critical field H c2 and 
on the thermal conductivity at very low temperatureiS. 
No "universal limit" has been observed, a surprising re- 
sult given the fact that the thermal conductivity of the 
pure sample is perfectly understood with the most pop- 
ular scenarios22*2i which predict a universal behaviour 
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at least in the basal plane. But measurements of the 
specific heat of irradiated samples already indicated a 
very unexpected behaviour, with a large broadening of 
the transition^. Theoretical investigations of the specific 
heat broadening due to the natural (gaussian) dispersion 
of defect concentrations 2 ^ cannot explain the importance 
of the effect observed, and further experimental charac- 
terisation of the induced disorder seemed necessary. 

In this paper, we present new results on the specific 
heat and transport measurements of UPt 3 samples 
damaged by high energy electrons. We show that the 
sample thickness has a dramatic influence on the results, 
explaining the previous failure of the measurements^. 
We give quantitative estimates of the distribution of 
defects in the sample, and determine the conditions 
for a clean test of the universal limit in this system. 
The experimental procedure and the main results are 
described in section II. In section III, an introduction 
to the creation of defects by electron irradiation is 
presented. Section IV presents a phenomenological 
model used to describe quantitatively our results. In the 
last sections, the previous results are briefly discussed in 
the framework of our description. 
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FIG. 1: Temperature dependence of the resistivity of the 
three samples of UPt3, Rl R2 and R3 (with a thickness of re- 
spectively 190, 90 and 70 /im) after an irradiation by electrons 
of energy of 2.5 MeV with a fluence of 3.5 10 18 xe"/cm 2 . The 
sample R2 is also shown before irradiation (closed circles) and 
is typical of pure samples of UPt3. The current is oriented 
along the a-axis except for sample Rl (c-axis). The left scale 
respects the ratio of the anisotropy of the effective mass of 
the quasi-particules observed in pure samples (2.5). 



II. EXPERIMENTAL DETAILS AND MAIN 
RESULTS 

Pure samples were cut from the same single crys- 
tal grown under ultra-high vacuum by the Czochralski 
method. The shape of the sample for resistivity mea- 
surements was a bar of typical length 5 mm and width 
0.23 mm. The samples Rl, R2 and R3 have been thinned 
down to 190, 90 and 70 /im respectively, in order to probe 
the depth of the irradiation damage. R4 and R5 with a 
respective thickness of 220 /im and 240 /im were prepared 
to study the evolution of the normal state transport prop- 
erties for different rates of irradiation. Rl and R4 were 
cut with the electrical current along the c-axis, R2, R3 
and R5 with the electrical current in the basal plane. The 
two samples measured by specific heat (CI and C2) were 
plates 270 /im thick and a surface roughly 2x2 mm 2 . Af- 
ter being cut and/or thinned, all samples were annealed 
under ultra high vacuum ( < 3.10~ 10 Torr) at 950°C dur- 
ing 5 days. The quality of the resistivity samples has not 
been altered by the thinning process since the three sam- 
ples exhibited a high critical temperature T c of 0.55 K 
(Rl) and 0.56 K (R2 and R3). 

The irradiations by high energy electrons have been 
performed in the low temperature facility of the Van 
der Graaff accelerator of the "Laboratoire des solidcs 
irradies" of Ecole Polytechnique at Palaiseau (France). 
The resistivity samples were held next to each other on 
a thin fibre glass plate with GE varnish. The thinnest 
dimension was parallel to the electron beam. This plate 
was wrapped in a thin copper foil, held rigidly inside a 
liquid hydrogen bath. The three samples were irradiated 



simultaneously, and they remained well within the ho- 
mogeneous region of the electron beam. Electrical con- 
tacts on the samples were done with silver paste (sol- 
der might have reannealed part of the defects), removed 
before, and made again after the irradiation. The sam- 
ples CI and C2 for the specific heat measurements were 
wrapped directly the thin copper foil, and were irradiated 
separately. Incident electrons were accelerated at an en- 
ergy of 2.5 MeV for samples for resistivity measurements 
and CI, and with an energy of 1.5 MeV in the case of 
sample C2. The irradiation takes place in liquid hydro- 
gen. Two Faraday cups located in the front and behind 
the samples arc used to estimate the irradiation fluence 
at the level of the sample^. The electron flux is limited 
to ~T0 14 e~/cm 2 /s in order to avoid the heating of the 
sample during irradiation. Fluences ($) between 1.8 and 
13.9 x 10 18 electron/cm 2 (e~/cm 2 ) have been used. After 
irradiation, the samples were warmed up to room temper- 
ature and transferred in different experimental set-ups to 
perform resistivity and specific heat measurements. This 
procedure results in some annealing of the defects created 
at low temperature. 

The resistivity measurements have been performed 
with a traditional four wires technique with a 17 Hz AC- 
current bridge. The specific heat measurements have 
been performed with a semi adiabatic technique. All 
samples have been measured before and after irradiation. 

When electrons go through the matter, a part of their 
energy is lost via inelastic collisions with the target elec- 
trons. Typically, the energy loss is of the order of 1-2 
keV//im2£. In the case of 100-200 /im thick the energy 
loss can be as high as 400 keV. But as long as the energy 
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transferred to the target nucleus, which depends on the 
atomic weight of the nucleus and on the incident energy 
of the electron, is much larger than the minimum energy 
required to displace it (threshold energy), the cross sec- 
tion for the atomic displacement does not depend drasti- 
cally on this energy. So, in this case, the decrease of the 
electron energy through the target compound will not in- 
duce strong inhomogeneities in the damage distribution. 
However, if the energy transferred to the target nucleus is 
close to the threshold energy, one could expect that the 
decrease of the incident energy along the path induces 
that electrons will not be longer able to displace atoms 
in the whole thickness of the samples. In this latter case 
the critical temperature can be expected to change dif- 
ferently along the electron path depth. We suspect that 
this situation could explain the large broadening of the 
specific heat transition observed previously on irradiated 
UPt3 samples^. In order to test this assumption, we 
have performed two kinds of measurements : (i) we have 
used lower energy electrons (if close to the threshold, the 
irradiation effects should then be drastically reduced), 
and (ii) we have measured irradiated samples of various 
thicknesses in order to check a possible dependence of T c 
with thickness. 

Figure ^ shows the result of the last test. It displays 
the resistivity of the three samples of different thickness, 
irradiated during the same experiment and with the same 
fluence <t> = 3.5 x 10 18 e~/cm 2 . The most striking result 
is that indeed, the thinner the sample, the stronger the 
drop in T c . This correlation between sample thickness 
and T c change is a direct proof of the inhomogeneity of 
the irradiation damage. As regard the normal state prop- 
erties, the temperature dependence of the resistivity re- 
mains quadratic (p(T) = po+A.T 2 ) for all samples. Both 
terms po and A increase under irradiation, and again, po 
increases more strongly when the sample is thinner. How- 
ever, A is closer to the value of a pure sample when the 
irradiated sample is thinner. We take this as pointing 
to an increased average disorder for thinner samples and 
a wider distribution of defects for thicker samples. We 
shall make these conclusions more quantitative in section 
V. 

Figure |2] presents the effects of irradiation at two dif- 
ferent energies on the specific heat. The top panel com- 
pares virgin sample C2, and the same sample irradiated 
with a fluence of 4.5 x 10 18 e~/cm 2 at 1.5 MeV : within 
experimental accuracy, there is no effect at all of the ir- 
radiation, both in the normal and superconducting state. 
The simplest interpretation is that with an incident en- 
ergy of 1.5 MeV, the transmitted energy to the uranium 
and platinum atoms is below the threshold energy for 
defect formation. These results at 1.5 MeV should be 
contrasted with the bottom panel of figure |21 It shows 
the drastic effects at 2.5 MeV2£ on sample CI, with a new 
measurement at a low fluence of 1.8 x 10 18 e~/cm 2 . The 
virgin sample CI (like C2) exhibits the well known double 
specific heat jump at temperatures (defined by an onset 
criteria) of T c +=550 mK and T c -=498 mK respectively. 
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FIG. 2: Temperature dependence of the specific heat of the 
same sample of UPt3 before (open circle) and after damaged 
by irradiations. The irradiations were performed with in- 
cident electrons at an energy of 1.5MeV (top panel sample 
C2) and 2.5 MeV (bottom panel sample CI). The top panel 
demonstrate that at 1.5MeV, no effects are observed on the 
specific heat (no defects creation). On the bottom panel, the 
line for the sample before irradiation is a phenomenological 
description of the specific heat behaviour with a small homo- 
geneous distribution of T c of 16mK. After irradiation, a good 
fit in the whole temperature range and for the two fluence is 
obtained with an exponentional decrease of the rate of defects 
in the depth of the sample (see text). The exponential dis- 
tribution is characterised by the length A = 90 pm and the 
decrease of T c for an infinitely thin sample by the parameter 
dT c /d$=-115 mK/10 18 e"/cm 2 . 



Both transitions show only a small broadening of 16 and 
12 mK respectively, demonstrating the very good quality 
of the sample. Note that at very low temperature the 
strong upturn in the specific heat already reported^ is 
well observed. The origin of this upturn could be a mag- 
netic ordering transition but is not directly connected to 
the superconducting state. In the following analysis this 
upturn will be ignored. After irradiation, the same sam- 
ple shows a strong broadening of the specific heat : the 
two superconducting transitions are already hardly dis- 
tinguishable for a fluence of 1.8 x 10 18 e~/cm 2 , and are 
totally smeared at a higher fluence. 

By contrast, the onset critical temperature shifts only 
very little, much less than observed by resistivity on 
thin samples (R3 for example) . From this point of view, 
the irradiation effects are different from the annealing 
effects where broadening and T c decrease are intimately 
correlated^. It points again to an inhomogeneous 
distribution of point defects. Let us note that normal 
state properties arc little affected by irradiation (apart 
from a smaller mean free path). Indeed, above the 
superconducting transition, the sample exhibits a clear 
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linear dependence of the specific heat (C = jT) with a 
large value of the Sommerfeld coefficient 7 characteristic 
of the heavy fermion Fermi liquid regime. The irra- 
diations do not affect this Fermi liquid behaviour nor 
the value of 7. They are also reversible : the original 
specific heat has been recovered after a subsequent an- 
nealing at 700° C during 5 days under ultra high vacuum. 



III. CREATION OF DEFECTS BY HIGH 
ENERGY ELECTRON IRRADIATION 

We have stated that the experimental behaviours of 
the resistivity and the specific heat suggest that 2.5 MeV 
is very close to the energy necessary to create stable de- 
fects in UPt3. In this section, we will present a short 
introduction to the creation of defects by high energy 
electron irradiation with a special emphasis to the spe- 
cial case of UPt3. A plausible qualitative explanation 
of the origin of such inhomogeneities is proposed. We 
then develop a simple model and introduce a character- 
istic length of defect creation which allows us to describe 
phenomenologically and quantitatively the broadening of 
the specific heat transition. This model is then applied 
to the resistivity results. 

The process of defect creation in a sample irradiated 
by high energy electrons is based on elastic collisions be- 
tween the incident electrons and the atomic nuclei of the 
crystal^. There is an activation energy - the threshold 
energy E t - necessary to create stable vacancy-interstitial 
defects (Frenkel pairs), i.e. vacancy and interstitial sepa- 
rated enough to avoid spontaneous recombination. In the 
case of monoatomic hexagonal compact structures, E t is 
typically about a few ten eV— . This value has to be com- 
pared to the maximum energy transferred to the atom by 
an incident electron which is inversely proportional to the 
mass of the nuclei. In UPt3, Uranium and Platinum are 
heavy nuclei with Z=238 and Z=I95 respectively. The 
maximum transmitted energy is respectively of 79.5 and 
97.0 eV for incident electrons of energy 2.5 MeV, and 34.1 
and 41.7 eV for an incident energy of 1.5 MeV. Therefore 
the fact that irradiation with 1.5 Mev electrons has no 
effect on the specific heat gives us a bottom limit for the 
threshold energies in UPt3 equal respectively to 34.1 and 
41.7 eV for U and Pt atoms. These values are quite high 
and might be explained by the high compactness of the 
crystallographic structure in UPt 3 . 

The cross section for displacement of an atom by an 
electron with an incident energy E depends both on 
the threshold energy and on the energy transferred to 
the atoms. For thick samples, it is then necessary to 
take into account the energy losses along the electron 
path. It is found to be equal to -2.37 keV/^m and to 
remain approximately constant over a depth of 270 /im 
(corresponding to the thickest sample of this study). 
This decrease of incident energy as a function of the 
penetration depth inside the sample induces a decrease 



of the cross section for elastic collisions and therefore 
of the concentration of defects. This is a schematic 
mechanism from which the rate of defects creation 
depends on the depth inside the sample. It is even 
stronger if the initial energy is barely enough for the 
transmitted energy to exceed the threshold. A detailed 
rigorous description of rate of defect creation inside the 
sample as a function of energy is certainly possible. Such 
analyses have been effectively performed, for example 
to determine the threshold energies in the diatomic 
compound TaC— 4. But it is not our purpose here, which 
is rather, first, to check that this is indeed responsible 
for the measured effects and, second, to determine which 
experimental conditions have to be used in order to 
prevent such an inhomogeneous distribution. This is the 
aim of the phenomenological model developed in the 
next section. 



IV. PHENOMENOLOGICAL MODEL 

In order to quantify our hypothesis, we need the con- 
centration of defects as a function of depth in the sample 
(c(x), x = corresponding to the side of the sample fac- 
ing directly the incident beam) . As explained above, c{x) 
which is of course proportional to fluence, is not known 
exactly and depends on parameters like the energy de- 
pendence of the cross section and the decrease of the en- 
ergy of the incident electrons due to inelastic scattering. 
So we have chosen to use a simple phenomenological ex- 
ponential dependence, parameterised by a characteristic 
length for defect creation A : c(x) = c(0)exp(— x/X). 

The second element is the relationship between defects 
concentration and T c decrease. This is well known and 
studied, both experimentally and theoretically. Indeed, 
it has been already demonstrated that the critical tem- 
perature of UPt3 is strongly dependent on the concen- 
tration of defects 13-1418 . This feature is common to all 
unconventional superconductors like high-T c oxides (see 
for example^) or Si^RuO^SS. The decrease of the critical 
temperature in a conventional or unconventional super- 
conductor due to a pair breaking mechanism is always 
described by the Abrikosov-Gor'kov formula^: 

T c is the critical temperature, T c q the critical tempera- 
ture of the sample without pair breaking, ip the digamma 
function and a characterises the strength of the pair 
breaking mechanism. In the limit of small pair break- 
ing (a — > 0), the decrease of the critical temperature is 
linear and equal to ■rr L - Larkin has shown that for an 
unconventional superconductor with isotropic impurity 
scattering, a is equal to 4p where V is the scattering 
rate^i. At fixed incident energy E, for an infinitely thin 
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sample, the scattering rate is proportional to the num- 
ber of defects and therefore to the fluence. As the cross 
section of individual scattering centers is not known, we 
rather define a new parameter : | (E) | , the initial rate 
of decrease of T c with fluence, for an infinitely thin sam- 
ple. For a given material, this parameter will only depend 
on the (fixed) incident energy. Then, the local critical 
temperature T c {x) at depth x in the sample is given by 
equation JQ) with : 
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A strong test of the model is to reproduce the specific 
heat behaviour. Again, exact calculations of the specific 
heat for a given homogeneous concentration of defects are 
possible^. They depend on the symmetry of the order 
parameter, but they do not account for the double tran- 
sition. So we used a much simpler approach, which ap- 
pears more proportionate to the ambition of our model. 
We first make a simple fit of the specific heat of the virgin 
sample, which defines a function Ct c0 (T) which respects 
the entropy conservation. For a homogeneous sample of 
reduced temperature T c (due to an homogeneous distri- 
bution of pair breaking centers), the specific heat will 
be Ct c (T). It has a residual linear term, due to impu- 
rity bound states^. This term is estimated crudely as 

= 7 Tc ° ~ Tc , where 7 is the Sommerfeld 



7r 



C(T) 
T 



coefficient in the normal state. We then simply rescale 
with respect to T c the specific heat of the virgin sample, 
using the following expression which respects entropy bal- 
ance : 
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A simple integration over the sample thickness using 
equations El and ^ gives the fit of the specific heat of 
irradiated samples for each fluence, which depends only 
of the two free parameters A and | (E) | . On figure [21 
the full lines are the results of the fit for A = 90/im and 
\^(2.5MeV)\ = 115mK/10 18 e-.cm- 2 . It is gratifying 
to see that it reproduces well both the broadening and 
the small shift of the onset critical temperature for both 
fluences and with the same parameters. More insight 
into the physical effects at work is gained looking at the 
distribution of T c inside the sample deduced from these 
parameters. 

Let us first discuss the T c change as measured by 
resistivity. As soon as a thin layer of the sample is 
superconducting, it will short circuit the sample. In 
the present geometry of the irradiation, the part of the 
sample with the largest critical temperature (producing 
the drop to the zero resistivity on figure is the side 
opposite to the electron beam. It should be the same 
part of the sample that gives the onset T c for the 



FIG. 3: Dependence of the critical temperature (T c ) normal- 
ized to the critical temperature before irradiation (T c o)with 
the thickness of a sample irradiated with electron a fluence 
$ = 3.5 10 18 e"/cm 2 . The closed circle are the experi- 
mental data measured by resistivity. The solid line is the 
value expected in function of the depth for an exponential 
dependence of the rate of defects with the parameters (see 
text) A = 80 jLtm and dT c /d$=-90 mK/10 18 e~/cm 2 . The 
dashed line is with the parameters A = 90 /im and dT c /d$=- 
115 mK/10 18 e~/cm 2 used for the specific heat fit on figure 2. 
The inset is a schematic description of the geometry of the 
irradiation of the samples. 



specific heat. On figure the ratio between the critical 
temperature after and before irradiation with the same 
fluence for samples of different thickness is reported. 
The solid curve is T c (a;), where x is the sample thickness, 
as deduced from equations ^ an( i 121 The parameters 
differ slightly than from those of the specific heat fit 
(A = 80^m and |^(2.5MeV)| = 90mK/10 18 e^.cm- 2 ). 
For comparison, the dashed line shows the expected 
T c variation for the resistivity experiment using the 
parameters of the fit of C' p . However, according to 
the crudeness of our model, and the difficulties always 
met in the comparison of T c between specific heat 
and resistivity, we believe that it is a "good enough" 
quantitative agreement. Of course, the lines in figure |3| 
also represent the distribution of T c as predicted by our 
model in a single sample. For the intermediate fluence 
used in the resistivity experiment, it can be seen that 
the distribution of T c is very large and extends very 
quickly from almost T c o (small change in the T c onset 
for a sample 3 times thicker than A) to 0. The amazing 
broadening of the specific heat transition is therefore 
simply explained by a trivial effect of sample geometry 
and small A, due to "unlucky" experimental conditions 
(2.5MeV being too close to the threshold for defect 
creation) . 



CONSEQUENCES ON THE TRANSPORT 
MEASUREMENTS 



Let us summarize some further consequences on trans- 
port measurements. The fact that the Sommerfeld coeffi- 
cient 7 of the normal state specific heat does not depend 
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on the irradiation shows that the density of states at the 
Fermi level is not affected. This effect seems to contradict 
the increase of the quadratic term of the temperature de- 
pendence of the resistivity (AT 2 ) under irradiation (see 
section I) . Indeed A and 7 are related by the well known 
Kadowaki- Woods ratio^. One can then speculate that 
the apparent increase of A is an artifact of the distri- 
bution of defects. Indeed, when the whole sample is in 
the normal state, the measured conductivity is an av- 
erage of the conductivity of the different layers of the 
sample, which are all in parallel owing to the contact ge- 
ometry. In the simplest models which assume that the 
Matthiessen's rule is verified, only the residual resistivity 
should dependent on the impurity concentration. Then, 
the temperature dependence of the conductivity should 
be described by : 

where t is the thickness of the sample, po(x) the 
residual resistivity at a depth x, and AT 2 the intrinsic 
quadratic term, independent of defects concentration. In 
general the calculation has to be done numerically, how- 
ever in the limit where AT 2 is small compared to po(x), 
the temperature dependence of the resistivity deduced 
from equation 0] will be of the form 1/ctq + AT 2 with : 
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FIG. 4: Dependence of the residual resistivity extrapolated 
to T=0K (bottom) and of the quadratic temperature depen- 
dence of the resistivity (top) for different rates of irradiation 
for the two samples R4 and R5 of respective crystallographic 
orientations : current along c-axis and current along a-axis. 
R4 (squares) is 220/im thick, and R5 (circles) is 240/im thick. 
Lines are the dependence calculated for a rate of defects de- 
creasing exponentially in the thickness of the sample with a 
characteristic length of QOfim (plain) or 80^im (dashed). Only 
the residual resistivity is assumed to be modified by the de- 
fects, with a rate of 1.4 (resp. 0.6) pfl. cm/10 18 xe~.cm -2 for 
the current along the a (resp. c) direction. 
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where (...) is the average over the depth. An increase of 
the dispersion of the defects concentration can therefore 
be at the origin of an apparent increase of the quadratic 
temperature dependence. Qualitatively, such an effect is 
observed with the increase of the thickness of the sample, 
but also for an increase of the fluence for the same sample. 

Quantitatively, since the dependence of the residual re- 
sistivity for an infinitely thin sample is proportional to 
the density of defects, again, a unique parameter ^£f~ 
(where i represents the two different crystallographic ori- 
entations) controls the increase of the apparent residual 
resistivity and of the quadratic term of a thick sample in 
each direction. At the temperature T, the resistivity at 
depth x inside the sample irradiated with a fluence $ is 
given by: 



Po(x) = P o + ^*exp(-x/\) + AT 2 (7) 

where po is the residual resistivity before irradiation 
and A the characteristic length already defined above. 
The result of this calculation with a parameter of 
1.4 (resp. 0.6) /xf2.cm/10 18 e~/cm 2 for a current along 



the a-axis (resp. c-axis) is shown on figure 4 for the 
same sample after different rate of irradiation and on 
the table 1 for a same fluence on samples of different 
thickness. The general trend is well described. 

Of course, the most interesting transport measurement 
is the thermal conductivity k in the superconducting 
state. The measurements of k on irradiated samples 



Thickness 


PO 
p Pure 


A 

A Pure 


(urn) 


meas. (Calc.) 


meas. (Calc.) 


70 


9.3 (12) 


1.12 (1.03) 


90 


8.7 (10.6) 


1.15 (1.05) 


190 


6.9 (5.6) 


1.30 (1.12) 



TABLE I: Ratio after and before irradiation by electrons of 
energy of 2.5 MeV with a fluence of 3.5 10 18 e~/cm 2 of (a) 
the residual resistivity extrapolated to T=0K and (b) the 
quadratic temperature dependence of the resistivity. The val- 
ues measured are compared with the values numerically cal- 
culated (between bracket) with an exponential dependence 
of the concentration of defects in the thickness of the sam- 
ples. We have also assumed that only the residual resistiv- 
ity changes (for an homogeneous distribution of defects, see 
text). The parameters used are X=80pm (see text) and a vari- 
ation of the resistivity 1.4 (0.6) p£l. cm/10 18 e~/cm 2 in the a 
(c) direction. These parameters are the same as those used 
in the figure 4. 
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reported earlier by our groupi* have been performed 
on samples typically of several hundreds \im thick. 
The irradiation parameters were an incident energy of 
2.5MeV and a fluence of 3.17x 10 18 e"/cm 2 and llxlO 18 
c~/cm 2 . The main result was the observation of a 
dramatic decrease of the thermal conductivity in the 
superconducting phase and a large increase of the zero 
temperature extrapolation of k/T with the fluence of 
the irradiation. The huge distribution of T c induced by 
both irradiations in these experimental conditions (down 
to T c = K) give a simple and natural explanation to 
these features. So the question of a "universal limit" in 
UPt 3 remains completely open, and should be settled by 
measurements on very thin samples irradiated on both 
sides to minimize the defects inhomogeneities. From 
the parameters used in our model, we estimate that 
for samples 70/zm thick, irradiated on both faces with 
total fluences below 1.5x 10 18 e~/cm 2 , the variation of 
T c should remain below 20% (from T c =0.55 K before 
irradiation to 0.46 K after), with an inhomogeneous 
broadening below 10% of the variation of T c (i.e. 10 mK 
which is small compare to the broadening of the "pure" 
samples). Let us note that for the same experiments 
performed on the bismuth compounds ' , the samples 
were already much thinner (« 20/j,mM?), the transmitted 
energies twice larger (due to lighter elements) and 
the threshold energy 4 times smaller—, ensuring that 
damage was homogenously distributed through the 
sample in contrast to the situation encountered for 
UPt3. The same "good" conditions should also be found 
in the compounds composed of lighter elements like 



Cerium based heavy fermion systems or borocarbide 
superconductors. 



VI. CONCLUSION 

In conclusion, we have measured and analyzed the 
specific heat of thick samples of UPt 3 damaged by 
high-energy electron irradiation. A strong smearing of 
the specific heat jump at the superconducting transition 
and a large increase of the residual C/T term are ob- 
served. We have shown that a strong dependence of the 
concentration of defects along the sample thickness can 
explain quantitatively this feature. Such inhomogeneities 
have a natural explanation with the heavy elements of 
UPt3. This is further confirmed by the dependence of 
the critical temperature measured by resistivity with 
the thickness of the sample. Such macroscopic inhomo- 
geneities give also a natural explanation to the increase 
of the quadratic term of the temperature dependence 
of the resistivity, and to the apparent absence of a 
universal limit in the heat transport measurements 
reported previously. 
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